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We discuss the influence of the cosmological constant on the gravitational equations of motion 
of bodies with arbitrary masses and eventually solve the two-body problem. Observational con- 
straints are derived from measurements of the periastron advance in stellar systems, in particular 
binary pulsars and the solar system. Up to now, Earth and Mars data give the best constraint, 
A < 10 -36 km -2 ; bounds from binary pulsars are potentially competitive with limits from inter- 
planetary measurements. If properly accounting for the gravito-magnetic effect, this upper limit on 
A could greatly improve in the near future thanks to new data from planned or already operating 
space-missions. 
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I. INTRODUCTION 



Aged nearly one century, Einstein's cosmological constant A still keeps unchanged its cool role to solve problems. 
A, despite being just one number, was able to respond to very different needs of the scientific community, from 
theoretical prejudices about universe being static (which provided the original motivation for introducing A in 1917) 
to observational hints that the universe is dominated by unclustered energy density exerting negative pressure, as 
\^ , required by data of exquisite quality which became available in the last couple of decades. Although it is apparently 
plagued by some theoretical problems about its size and the coincidence that just in the current phase of the universe 
the energy contribution from A is of the same order of that from non-relativistic matter, the cosmological constant still 
provides the most economical and simplest explanation for all the cosmological observations [l|. The interpretation 
of the cosmological constant is a very fascinating and traditional topic. A might be connected to the vacuum density, 
as suggested by various authors (see for an historical account), and could offer the greatest contribution from 
! ; cosmology to fundamental physics. 

The big interest in the cosmological constant has recursively raised attempts in putting observational bounds on its 
absolute value from completely different phenomena. A, supposed to be ~ 10~ 46 km~ 2 from observational cosmology 
analyses, is obviously of relevance on cosmological scales but it could play some role also in local problems. Up to 
C$ • now, no convincing methods for constraining A in a laboratory have been proposed [3j, but interesting results have 
been obtained considering planetary motions in the solar system |4j, y, |6j . The effects of A become stronger for diluted 
mass conglomeration but they get enhanced also through various mechanisms 0, • As an example, conditions for 
' the virial equilibrium can be affected by A for highly flattened objects 0. On the scale of the Local Volume, a 
cosmological constant could have observable consequences by producing lower velocity dispersion around the Hubble 
flowg. 

Up to now, local physical consequences of the existence of a cosmological constant were investigated studying the 
motion of test bodies in the gravitational field of a very large mass. This one-body problem can be properly considered 
in the framework of the spherically symmetric Schwarzschild vacuum solution with a cosmological constant, also known 
as Schwarzschild-de Sitter or Kottler space-time. The rotation of the central source can also be accounted for using 
the so-called Kerr-de Sitter space-time Here, we carry out an analysis of the gravitational A^-body problem with 
arbitrary masses in the weak field limit with a cosmological constant. This study is motivated by the more and more 
central role of binary pulsars, from the discovery of the pulsar PSR B1913+16 in 1974 [l(J], in testing gravitational 
and relativistic effects. The gravitational two-body equations of motion for arbitrary masses were first derived in 
absence of spin by Einstein, Infeld and Hoffmann (EIH) [ll|. The problem was later addressed in more general cases, 
subsequently accounting for spins and quadrupole moments |l2t and reference therein]. Here, we take the further step 
to consider a cosmological constant. 
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The paper is as follows. In section [D] we discuss the gravitational weak field limit in presence of a cosmological 
constant and introduce the relevant approximations. Section 11111 presents the generalization of the EIH equations 
of motion, whereas section IIVI is devoted to the study of the two-body problem. In section ^ we review how 
measurements of precession of pericentre in stellar system can constrain A. In particular, we consider binary pulsars 
and the solar system. Section IVTl contains some final considerations. 



II. FIELD EQUATIONS WITH A COSMOLOGICAL CONSTANT IN POST-NEWTONIAN 

APPROXIMATION 



Einstein's equations with the cosmological constant are 



PA— 87rG . . 

R/iu — J^9fj.u — ^ 4 &fj,v, (1) 

where G is the gravitational constant, c the vacuum speed of light and 

S^u — T^u — —g^uTx, (2) 

with being the energy- momentum tensor. The weak field expansion can start by introducing a nearly Lorentzian 
system for weak, quasi-stationary fields, in which 

9nv = Vfiu + V" |^| < 1- (3) 

Actually, the Minkowski metric rj^v is not a vacuum solution of the field equations with a cosmological constant, 
but for |A| <C 1 an approximate solution in a finite region can still be found by an expansion around rj^ v . In the 
post-Newtonian (pN) approximation, metric components can be expanded in powers of 

c J9i£)'\^i, (4) 

\ &R J c p 

where M, R, v, p and p represent typical values for the mass, length, velocity, pressure and energy density of the 
system, respectively. In what follows, ^ n 'g^ u and will denote terms of order e™ and e n (M/R), respectively. 

To perform a proper treatment in presence of a A-term we have to consider the suitable approximation order for 
A. We assume that the size of the contributions due to the cosmological constant is at most comparable to the 
post-Newtonian terms, i.e, O(Agoo) ^ 0(G ( 2 )T 00 /c 2 ). This condition can be rewritten as 

A < §, (5) 

where R g = GM/c 2 is the gravitational radius. Eq. © is easily satisfied by gravitational bound systems with M ~ M Q 
and R ~ 1 — 10 2 AU if A < 10~ 33 km, a value well above the estimated one from cosmological constraints and also 
greater than the limits we will derive in section Ivl considering stellar systems. Hereafter, we will put c = 1. With 
such an approximation order, we can use classical results within the standard pN gauge. Following Straumann |l3j |. 
the approximate field equations read 

A( 2 >. 900 = -SttG ^T 00 , (6) 

A (4) 5oo = W 9ij (2} 9oo^ + (2) gvj {2) 9oo, ~ \ {2) 9oo, (2) 5oo, - \ {2) 9oo, < a W (7) 

- 8itG ( ^T 00 - 2 ( 2 > 500 (0) T 00 + ^ 1*) + 2A, 

A^g ot = -\ ^g jjfii + (2) 9ijfij + !6nG ^T i0 , (8) 
A^ 9ll = -SnGSij ^T Q0 . (9) 



The components of the metric can be expressed in terms of potentials. Let </>n be the Newtonian potential, 

r (o)r 00 (t x'l 

<f> N = -G / , l ', (10) 
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According to our approximation order, the cosmological constant appears only in the equation for ^'goo. This can be 
re-arranged to give 

AW (300 + 2$,) = -8nG ( (2) T 00 + (2) T**) + 2A (11) 
Together with the classical pN potential i/j, 

^ = -G f r^-^j ( (2) T 00 , (12) 

./ x - x V / 



we introduce <f>A, solution of 



A0 A = -A. (13) 



In presence of a cosmological constant, there is an upper limit on the maximum distance within which the Newtonian 
limit holds and boundary conditions must then be chosen at a finite range . When these boundary conditions are 
chosen on a sphere whose origin coincides with the origin of the coordinate system, can be expressed as 

A = -pA|x| 2 , (14) 
6 

where we have neglected correction terms which appear because of boundary conditions. Due to a positive cosmological 
constant, the origin of the coordinate system has a distinguished dynamical role with a radial force directed away 
from it |l5j . Since the choice of the origin is arbitrary, any point in the space will experience repulsion from any other 
point. Finally, introducing the standard pN potentials, 

6 = - 4G /u^m (1) r i0 (t,x'), (is) 

X = -f J |x-x'|<°>r°°dV(t,x'), (16) 

the metric components read 

Wgoo = -2<fa, (17) 

Wg 00 = -2(^ + V + ^a), (18) 

(2) 5i , = -26 ij( j> N , (19) 

(3) 5o* = Zi + X,io- (20) 

For a point-like mass at the centre of the coordinate system, the above expressions reduce to the weak field limit at 
large radii of the Kottler space-time. 

A. Equations of motion for a test particle 

The motion of a particle in an external gravitational field can be described by the Lagrangian 



V dt J \ dt 



Using the metric components in equations I|17l20fl . we get 



1 .2 ,1.4 j. 1 j.2 „,. j. 3 a ..2 , ..% ( i- , d X 



The corresponding Euler-Lagrange equations of motion in a 3-dimensional notation read, 

^ ~ -V (0 N + 2^ + V) - ^ - J^Vx + vx (VxO + 3v^ + 4v(v V)0 N - v 2 V^ N + ~x. (23) 
The above expression reduces to equation (20) in [Jj when neglecting pN corrections. 
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III. THE EINSTEIN- INFELD-HOFFMANN EQUATIONS 

Since the contribution from the cosmological constant is of higher-order, it does not couple with other corrections. 
The Lagrangian of an JV-body system of point-like particles can be written as 

C ~ £ (A=0) + 5C A , (24) 

where £a=o is the total Lagrangian in absence of A. The Lagrangian C a of particle a in the field of other particles is 

C a ~£ a(A =o) + -**, (25) 
where C a (A=o) 1S given in equation (5.94) in [l3T ] The total Lagrangian reads 

C ~ £(A= ) + |- m « X a' ( 26 ) 

a 

with £(a=o) given in ?13l . equation (5.95)]. The corresponding Euler-Lagrange equations are the Einstein-Infcld- 
Hoffmann equations corrected for a A term, 

*« = -° E (—) + ^ F PN(A=0) + ^x Q (27) 

Ma V 1 

where F p n(a=o) is the post-Newtonian perturbing function |l3l equation (5.96)]. 



IV. THE TWO-BODY PROBLEM 

The total Lagrangian for two particles can be written as 

C ~ -m & vl + G m& ™ h + -mbvl + 5C pN{A=Q) + 5C A (28) 

where x = x a — Xb is the separation vector and <5£ p n(a=o) an d 5C A are the pN and A-contributions, respectively. It 
is 

&C P N(A=o) = g(m a w a + m b u b ) + G— — [3(w a + u b ) - 7v a v b - (v a n)(v b n)J - — — (29) 

with n = x/x and 

SCa = ^{rn & x\ + m h xl). (30) 

Due to cosmological constant, the energy of the system is modified by a contribution —5C\. The pN and A corrections 
are additive and can be treated separately. We are interested in examining the effect of a non vanishing A term. Let 
us consider the centre of mass and relative motions. Introducing X = (m a x a + TO b Xb) /M , with M = m a + m^, the 
Lagrangian can be re-written as 

1 , A , 1 , A , Mu , s 

C ~ g + 6 MX + 2 M + 6 M + x' (31) 

with /i = m^m^/M . Due to cosmological constant, the centre of mass of the system is subject to an effective repulsive 
force given by AX/3 per unit mass. 

The equations for the relative motion are those of a test particle in a Schwarzschild-de Sitter space-time with a 
source mass equal to the total mass of the two-body system. Since the perturbation due to A is radial, the orbital 
angular momentum is conserved and the orbit is planar. The main effect of A on the orbital motion is a precession of 
the pericentre 0, ITfil and references therein] . Following the analysis of the Rung-Lenz vector in Q and restoring the 
c factors, we get for the contribution to the precession angular velocity due to A, 

^a = ^VT^, (32) 
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TABLE I: Limits on the cosmological constant due to extra-precession of the inner planets of the solar system. 



Name 


8uj a (arcsec/year) 


lua (deg/year) 


A Um (km" 2 ) 


Mercury 


-0.36(50) x 10~ 4 


9.61xl0 25 A/(l km -2 ) 


4xl0" 35 


Venus 


0.53(30) x 10~ 2 


2.51xl0 26 A/(l km" 2 ) 


9xl0 -33 


Earth 


-0.2(4) x 10~ 5 


4.08xl0 26 A/(l km" 2 ) 


lxl0~ 36 


Mars 


0.1(5) x 10" 5 


7.64xl0 26 A/(l km" 2 ) 


2xl0" 36 



"Prom 0. 



where e is the eccentricity and P\> the Keplerian period of the unperturbed orbit. This contribution should be 
considered together with the post-Newtonian periastron advance, i2> p n = 3(27r/P b ) 5 / 3 (GM/c 3 ) 2 / 3 (l - e 2 )" 1 . The 
ratio between these two contributions can be written as, 

w pN R g p 6i?| 

where R = a(l — e 2 ) 3 / 8 is a typical orbital radius with a the semi-major radius of the unperturbed orbit, p = 
M/(4irR 3 /3) is a typical density of the system and pa = c 2 A/8irG is the energy density associated to the cosmological 
constant. The effect of A can be significant for very wide systems with a very small mass. 



V. OBSERVATIONAL CONSTRAINTS 



In this section, we derive observational limits on A from orbital precession shifts in stellar systems and in the solar 
system. 



A. Interplanetary measures 

Precessions of the perihelia of the solar system planets have provided the most sensitive local tests for a cosmological 
constant so far 13, ■ Estimates of the anomalous perihelion advance were recently determined for Mercury, Earth 
and Mars p"?l Il8j. Such ephemerides were constructed integrating the equation of motion for all planets, the Sun, 
the Moon and largest asteroids and including rotations of the Earth and of the Moon, perturbations from the solar 
quadrupole mass moment and asteroid ring in the ecliptic plane. Extra-corrections to the known general relativistic 
predictions can be interpreted in terms of a cosmological constant effect. We considered the 1-er upper bounds. Results 
are listed in TableQ] Best constraints come from Earth and Mars observations, with A < 10 _36 km~ 2 . Major sources 
of systematic errors come from uncertainties about solar oblateness and from the grayito-magnetic contribution to 
secular advance of perihelion but their effect could be in principle accounted for 19] . In particular, the general 
relativistic Lense-Thirring secular precession of perihelia is compatible with the determined extra-precessions |19| . 
The accuracy in determining the planetary orbital motions will further improve with data from space-missions like 
BcpiColombo, Messenger and Venus express. By considering a post-Newtonian dynamics inclusive of gravito-magnetic 
terms, the resulting residual extra-precessions should be reduced by several orders of magnitude, greatly improving 
the upper bound on A. 

The orbital motion of laser-ranged satellites around the Earth has been also considered to confirm general relativistic 
predictions. Observations of the rates of change of the nodal longitude of the LAGEOS satellites allowed to probe 
the Lense-Thirring effect with an accuracy of ~ 10%, i.e. about half a milliarcsecond per year [2(|. Other proposed 
missions, such as the LARES/WEBER-SAT satellite [2l|. should further increase this experimental precision. In 
general, since effects of A become significant only for very dilute systems, even very accurate measurements of orbital 
elements of Earth's satellites can not help in constraining the cosmological constant. For a satellite with a typical 
orbital semi-major axis of about 12,000 km, in order to get a bound on A as accurate as those inferred from Earth 
and Mars perihelion shifts (i.e. A < 10~ 36 km -2 ), changes in orbital elements should be measured with a today 
unattainable precision of a few tens of picoseconds of arc per year, about six order of magnitude better than today 
accuracy. 
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TABLE II: Binary pulsars with known post-Keplerian parameter ill and corresponding limits on the cosmological constant. 
The identification of the companion is often uncertain. We refer to the original papers for a complete discussion. 



PSR Name 


P b (days) 


e 


Co (deg/year) 


diA (deg/year) 






A Um (km 2 ) 


ref. 


Double neutron star binaries 


J1518+4904 


8.634000485 


0.2494849 


0.0111(2) 


9.335xl0 24 A/ 


1 km" 


2\ 


2x10" 


29 


[22] 


B1534+12 


0.2736767 


0.420737299153 


1.755805(3) 


2.772xl0 23 A/ 


1 km" 


2 x 


1x10" 


29 


[23] 


B1913+16 


0.323 


0.617 


4.226595(5) 


2.838xl0 23 A/ 


1 km" 


2\ 


2x10" 


29 


[23] 


J1756-2251 


0.319633898 


0.180567 


2.585(2) 


3.510xl0 23 A/ 


1 km" 


2\ 


6x10" 


27 


[24] 


J1811-1736 


18.779168 


0.82802 


0.009(2) 


1.176xl0 25 A/ 


1 km" 


2\ 


2x10" 


28 


[25] 


J1829+2456 


1.176028 


0.13914 


0.28(1) 


1.300xl0 24 A/ 


1 km" 


2\ 


8x10" 


27 


[26] 


B2127+11C 


0.68141 


0.335282052 


4.457(12) 


7.168xl0 23 A/ 


1 km" 


2\ 


2x10" 


26 


[23] 


B2303+46 


12.34 


0.65837 


0.01019(13) 


1.037xl0 25 A/ 


1 km" 


2\ 


1x10" 


29 


[2JJ 


Neutron star/white dwarf binaries 


J0621+1002 


8.3186813 


0.00245744 


0.0116(8) 


9.288xl0 24 A/ 


1 km" 


2\ 


9x10" 


29 


[28] 


J1141-6545 


0.171876 


0.1976509587 


5.3084(9) 


1.881xl0 23 A/ 


1 km" 


2\ 


5x10" 


27 


[23] 


J1713+0747 


67.82512987 


0.0000749406 


0.0006(4)" 


7.573xl0 25 A/ 


1 km" 


2 x 


8x10" 


30 


[29] 


B1802-07 


2.617 


0.212 


0.0578(16) 


2.856xl0 24 A/ 


1 km" 


2\ 


6x10" 


28 


[27] 


J1906+0746 


0.085303(2) 


0.165993045(8) 


7.57(3) 


9.392xl0 22 A/ 


1 km" 


2\ 


3x10" 


25 


[30] 


Double pulsars 


J0737-3039 


0.087779 


0.102251563 


16.90(1) 


9.750xl0 22 A/ 


1 km" 


2\ 


1x10" 


25 


[23] 


Unknown companion 


B 1820- 11 


357.7622(3) 


0.79462(1) 


o.or 


2.425xl0 26 A/ 


1 km" 


2\ 


4x10" 


29 


[31] 



"Upper limit 



B. Binary pulsars 

Binary pulsars have been providing unique possibilities of probing gravitational theories. Relativistic corrections 
to the binary equations of motion can be parameterized in terms of post-Keplerian parameters |32|. As seen be- 
fore, the advance of periastron of the orbit, u>, depends on the total mass of the system and on the cosmological 
constant. In principle, because Keplerian orbital parameters such as the eccentricity e and the orbital period Pb 
can be separately measured, the measurement of Cj together with any two other post-Keplerian parameters would 
provide three constraints on the two unknown masses and on the cosmological constant. As a matter of fact for real 
systems, the effect of A is much smaller than tL> P N, so that only upper bounds on the cosmological constant can be 
obtained by considering the uncertainty on the observed periastron shift. We considered binary systems with mea- 
sured periastron shift, see Tabled The effect of A is maximum for B1820-11 and PSR J1713+0747. Despite of the 
low accuracy in the measurement of Co, PSR J1713+0747 provides the best constraint on the cosmological constant, 
A < 8xl0 _30 km~ 2 . Uncertainties as low as Suj > 10~ 6 have been achieved for very well observed systems, such as 
B1913+16 and B1534+12. Such an accuracy for B1820-11 would allow to push the bound on A down to 10 _33 km~ 2 . 

Better constraints could be obtained by determining post-Keplerian parameters in very wide binary pulsars. We 
examined systems with known period and eccentricity as reported in [3j| . The binary pulsar having the most favourable 
orbital properties for better constraining A is the low eccentricity B0820+02, located in the Galactic disk, with 
w A ~ 1.4xl0 27 A/(l knT~ 2 ) deg /days. For binary pulsars J0407+1607, B1259-63, J1638-4715 and J2016+1948, the 
advance of periastron due to the cosmological constant is between 7 and 9xl0 26 A/(l km -2 ) deg /days. All of these 
shifts are of similar value or better than the Mars one. A determination of uj for B0820+02 with the accuracy obtained 
for B1913+16, i.e. 5di > 1CT 6 deg /days would allow to push the upper bound down to 10~ 34 — 10~ 33 km~ 2 . 

VI. CONCLUSIONS 

We considered the TV-body equations of motion in presence of a cosmological constant. The impact of A on the 
two-body system was explicitly derived. Due to the anti-gravity effect of the cosmological constant, the barycentre of 
the system drifts away. The relative motion is like that of a test particle in a Schwarzschild-de Sitter space-time with 
a source mass equal to the total mass of the two-body system. The main effect of A is the precession of the pericentre 
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on the orbital motion. 

We determined observational limits on the cosmological constant from measured periastron shifts. With respect 
to previous similar analyses performed in the past on solar system planets, our estimate was based on a recent 
determination of the planetary ephemerides properly accounting for the quadrupole moment of the Sun and for major 
asteroids. The best constraint comes from Mars and Earth, A < 1 — 2 x 10~ 36 km -2 . 

Due to the experimental accuracy, observational limits on A from binary pulsars are still not competitive with 
results from interplanetary measurements in the solar system. Accurate pericentre advance measurements in wide 
systems with orbital periods > 600 days could give an upper bound of A < 10~ 34 — 10~ 33 km~ 2 , if determined with the 
accuracy performed for B1913+16, i.e. Soj > 10~ 6 deg /years. For some binary pulsars, observations with an accuracy 
comparable to that achieved in the solar system could allow to get an upper limit on A as precise as one obtains from 
Mars data. 

The bound on A from Earth or Mars perihelion shift is nearly ~ 10 10 times weaker than the determination from 
observational cosmology, A ~ 10~ 46 km -2 , but it still gets some relevance. The cosmological constant might be the 
non perturbative trace of some quantum gravity aspect in the low energy limit fij • A is usually related to the vacuum 
energy density, whose properties depends on the scale at which it is probed [if. So that, in our opinion, it is still 
interesting to probe A on a scale of order of astronomical unit. Measurements of periastron shift should be much 
better in the next years. New data from space-missions should get a very high accuracy and might probe spin effects 
on the orbital motion 0, . A proper consideration of the gravito-magnetic effect in these analyses plays a central 
role to improve the limit on A by several orders of magnitude. 

Note added. After submission of this work, L. Iorio |35j presented an analysis of solar system data similar to our 
results in section IV Al 
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